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A meaningful probability distribution for measurements of a quantum stress tensor operator can 
only be obtained if the operator is averaged in time or in spacetime. This averaging can be regarded 
as a description of the measurement process. Realistic measurements can be expected to begin 
and end at finite times, which means that they are described by functions with compact support, 
which we will also take to be smooth. Here we study the probability distributions for stress tensor 
operators averaged with such functions of time, in the vacuum state of a massless free held. Our 
primary aim is to understand the asymptotic form of the distribution which describes the probability 
of large vacuum huctuations. Our approach involves asymptotic estimates for the high moments of 
the distribution. These estimates in turn may be used to obtain estimates for the asymptotic form 
of the probability distribution. Our results show that averaging over a hnite interval results in a 
probability distribution which falls more slowly than for the case of Lorentzian averaging, and both 
fall more slowly than exponentially. This indicates that vacuum huctuations effects can dominate 
over thermal huctuations in some circumstances. 

PACS numbers: 03.70.+k,02.50.Cw,04.62.+v,05.40.-a 


I. INTRODUCTION 


It is well known that the vacuum fluctuations of a linear quantum field operator are associated with a Gaussian 
probability distribution. More precisely, it is an averaged field operator that has a well dehned probability distribution. 
Let x) be a field operator or a derivative of a field operator. The averaging can be conducted along a timelike 
curve (or other timelike submanifold), over a spatial volume, or over a spacetime volume (see, e.g. Ref. [1] for a 
discussion and references). Here we illustrate the case of temporal averaging with a real-valued sampling function 
frit), with characteristic width r which satisfies 


fr{t)dt = 1 . 


( 1 ) 


Define the averaged field operator by 

/ OO 

fr{t)dt. ( 2 ) 

-OO 

This operator has finite moments in the vacuum state, fin = which are those of Gaussian distribution; this may 
be summarized in the formula 


E 

n—O 


{ixy 


n: 


fin = (e*^"^) = e 2 




( 3 ) 


Thus the odd moments vanish, and the even moments grow as n\ By the Hamburger moment theorem (See, for 
example, Ref. m), these moments uniquely define the probability distribution to be a Gaussian 


P{ip) = 


exp 


‘fl 

'2a 


( 4 ) 
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Here the variance a = {ip^) is determined by both the functional form and width of frit)'- crucially, cr is finite provided 
the Fourier transform /(w) decays sufficiently fast as a; —>■ ood Some physical effects of time averaged quantum electric 
field fluctuations have been discussed in Refs. Em, in examples where the physical situation can define the sampling 
function. 

The probability distributions for quadratic operators, such as the energy density or other components of the stress 
tensor, were discussed in Refs. ElE]- Here it is an average of a normal ordered operator which has finite moments. In 
the case of the stress tensor for conformal field theory in two spacetime dimensions with a Gaussian sampling function 
in either space or time, the probability distribution is a shifted gamma distribution E] with the lower bound of its 
support given by the quantum inequality bound [7]. This includes the case of the energy density of a massless scalar 
field, for which the quantum inequality bound was found by Flanagan [S]. 

The corresponding problem for quadratic operators in four spacetime dimensions is more complicated. Here the 
sampling must be in time or spacetime, not space alone. ^ The case of Lorentzian sampling in time of the stress tensor 
was treated in Ref. E], where several features of the probability distribution were inferred from the calculation of a 
finite set of moments. (See also Ref. [13], in which some combinatorial results used in Ref. jH] were proved.) In this 
case, the moments grow as (3n)!, and neither the Hamburger moment condition nor the related Stieltjes condition is 
satisfied, so the moments are not guaranteed to define a unique probability distribution (and the available evidence 
suggests that they do not). However, the asymptotic behavior of the probability distribution for large arguments can 
be inferred from the moments. This allows us to estimate the probability of especially large fluctuations. For example, 
let u be the Lorentzian average on a time scale r of the electromagnetic energy density. In units where h = c = I, 
which will be used throughout this paper, the dimensionless measure of the averaged energy density can be taken to 
be a:; = ut'^. The asymptotic form of the probability distribution for large x is of the form 

P{x) ^ x~^ exp{—ax^^^) , (5) 


where a is a numerical constant of order one. The x^/^ dependence in the exponential indicates that large energy 
density fluctuations are more likely than one might have expected. Large thermal fluctuations are suppressed by an 
exponential in energy, and hence can be dominated at larger energies by vacuum fluctuations. 

Quantum stress tensor fluctuations can have a variety of physical effects, and have been discussed by several 
authors [I4H26] . These effects include radiation pressure fluctuations [14] , and passive quantum fluctuations of gravity, 
where the gravitational field fluctuates in response to matter stress tensor fluctuations. Much of the previous work 
on this topic has involved calculating mean squared effects using integrals of a stress tensor correlation function. An 
exception is a study of black hole and Boltzmann brain nucleation based on Eq. ([^ in Ref. E] . A deeper understanding 
of the probability distribution will allow further investigation of the effects of large fluctuations. 

The purpose of the present paper is to explore the probability distributions associated with sampling functions with 
compact support in time. These are functions which are strictly zero outside of a finite interval, but can nonetheless 
be infinitely differentiable. Such functions seem to be the most appropriate descriptions of measurements which begin 
and end at finite times, and avoid the long temporal tails of functions like the Lorentzian. Some examples of such 
functions will be treated in Sec.jTTJ and in particular it will be shown constructively that, for any 0 < a < I, there are 
real-valued and nonnegative smooth compactly supported functions whose Fourier transforms are of form exp(—/3|a;|“) 
up to a fractional error of order where /? > 0. The parameter a controls the behavior of the test function near 

the switch-on and switch-off regions; we also give a physical motivation for the use of such functions. The asymptotic 
behavior of the moments of quadratic operators sampled with these functions will be treated in Sec. |HI| and shown 
to have a dominant contribution that grows like r((np + 2)/a), where p depends on the operator in question (p = 3 
for the energy density). The resulting asymptotic form for high moments will be used in Sec. IV to infer the tail of 
the probability distribution for large arguments, and hence to estimate the probability for large fluctuations. The tail 


^ The same is true for averaging over spatial and spacetime volumes, but is not true for averaging over spacelike curves. 

^ See again Ref. more recent developments using microlocal analysis permit a more geometric framework that can be used to 
understand the reason spatial averaging is excluded. The basic issue is the ability to restrict distributions to submanifolds; such 
restrictions certainly exist if the wavefront set of the distribution contains no covector that can annihilate all the tangent vectors to the 
submanifold (see, e.g., the discussion of timelike curves in llOl l. In the case of a (derivative of) a free field, the relevant wavefront set 
contains only null covectors, which cannot annihilate the tangent vectors to a timelike curve or a spacelike hypersurface of codimension 
1. However, in the case of quadratic operators the relevant wavefront set contains not only null, but also timelike covectors, which can 
annihilate all the tangent vectors to a spacelike submanifold, so the restriction cannot be achieved by these means. There is a more 
heuristic argument for the inadequacy of spatial averaging alone. A plane wave expansion of a quadratic operator will contain terms 
proportional to e*(*'+*^ Modes for which k = — k' will yield a divergent contribution to the moments which is not suppressed by 
spatial averaging. This effect is the reason why there are no quantum inequalities for spatially averaged stress tensors in four spacetime 
dimensions m- Similar problems affect averages along null curves m- 
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takes the form P{x) ~ cpa;^ exp(— oa;“/^) for certain parameters cq, a and b, and it will be shown that the probability 
of large fluctuations can be much greater than that for the case of Lorentzian averaged operators. Our results and 
their possible applications will be discussed in Sec. |Vj Three appendices contain technical details on various estimates 
used in the text. 


II. COMPACTLY SUPPORTED SAMPLING FUNCTIONS 
A. General results and constructions 

The sampling functions studied in this paper will be smooth and compactly supported nonnegative functions, 
which model a sampling period of finite duration. At the outset it is worth recalling some facts concerning the Fourier 
transforms of such functions. Let / be smooth and compactly supported. Then its transform 

/ OO 

( 6 ) 

-OO 


is clearly bounded, |/(w)| < J^^dt\f{t)\, and as this conclusion applies equally to the derivatives of /, we deduce 
that, for any n e N, 


I/HI = 


1 + ~ 1 + a;2 


(7) 


for some constant C depending on / and n. This is the well-known statement that the transform decays faster than 
any inverse power. On the other hand, the transform cannot decay too rapidly: if there are positive constants A and 
a such that |/(w)| < Ae““l“l for all w, then the inversion formula 


/W = 


( 8 ) 


provides an analytic extension of / to the strip |3z| < a of the complex plane. As / vanishes with all its derivatives 
outside its support, it must therefore vanish identically by uniqueness of analytic continuation. However, Ingham [27] 
has shown that there are functions of compact support with transforms obeying /(w) = 0(exp(—|a;|? 7 (|a;|))) as |a;| —)■ oo 
for any positive function r\ decreasing monotonically to zero as a; —>■ oo and such that 


auj - < (X). 


LU 


(9) 


In particular, Ingham’s result applies to any function r]{uj) = /3a;“ ^ for 0 < a < I and /3 > 0 establishing the 
existence of smooth compactly supported functions / with /(w) = At this stage, we remark that our 

usage of the ‘big-0’ notation will coincide with that in Ref. [3S]: to write that g{s) = 0{h{s)) on some set S means 
that there exists a constant A > 0 such that |g(s)| < A|/i(s)| for all s S S' - in fact we will only ever use this notation 
in situations where h{s) is nonnegative, so one has |g(s)| < Ah{s). To write that ^g{s) = 0{h{s)) as s —>■ oo’ indicates 
that there exists sq > 0 such that g{s) = 0{h{s)) on [sq, oo). 

Ingham’s construction defines the desired test function / in terms of its Fourier transform, which is a carefully 
arranged infinite product of sine functions. The precise behavior of f{t) (which can be regarded as an infinite 
convolution of top hat functions) is not investigated. Instead, we give a different construction that is slightly more 
explicit. To start, let (p G be a smooth nonnegative and integrable function so that (p(t) = 0 for t < 0 and 

(p{t) > 0 for t > 0, with (p{t) —>■ 0 as t —>■ -boo. Choosing 5 > 0, we may construct a smooth function with support 
[—(5, (5] by setting 


H{t) = (p{t + S)<p{6 — t). (10) 

The choice of 6 can be used to tune the shape of H. For instance, if p has its first local maximum at 6, then elementary 
differentiation shows easily that H has a local maximum at t = 0. On the other hand, we can achieve a broader peak 
by taking S slightly beyond the maximum of p, at the expense of making t = 0 into a local minimum of H; this is 
what we will do in the specific examples below. 
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As ip is smooth and integrable, it has a Fourier transform, which decays rapidly. Hence the convolution theorem 
applies, giving 


ff (w) = e 


i5ijj 


dw' 

27r 


ip{ui — 




( 11 ) 


where we have used elementary results on Fourier transforms of shifts and reflections, i.e., J^[p{t + d)](a;) = e*“'^(^(a;) 
and J^[p{5 — t)](w) = w). For w > 0, it is convenient to rewrite this expression as 


H{uj) = 23? , /(w) = 


!-uj/2 


duj' 

1'K 


(p{uj + uj')(f{uj')e 


/\^2iSuj' 


( 12 ) 


which may be proved by splitting the integration region in Eq. 0 into (—oo,w/2) and (a;/2, oo) and making the 
changes of variable uj' —> —u)' and oj' ^ to' — uj in the resulting integrals, thus yielding 


r^iSu) 


f-u;l2 


rh / 


and e 


— i5uj 


dw' 


/-w/2 


— p{-i^')‘p{-uj - uj')e 


f\^—2i5ijj' 


(13) 


which are seen to be mutually complex conjugate on recalling that ip{—uj) = ^*{ 10 ) because p is real-valued. As H 
is real and even, the same is true of H, and therefore Eq. (12) can be modified to hold for all uj by replacing all 
occurrences of w on the right-hand side by |a;|. 

Our interest is in the asymptotic behavior of H as |w| —>■ 00 . The specific examples chosen for p below have the 
property that the derivative p' decays more rapidly than p, which itself decays rapidly. Under these circumstances, 
the factors of p{uj + io') in the definition of I{oj) may be taken as approximately constant over the effective region of 
integration, leading to the approximation 


H{uj) ~ 23? ( e^^^p{uj) r 


/-u;/2 


27r 


(14) 


(We use ~ for an uncontrolled or heuristic approximation, reserving the symbol ^ for a rigorously established asymp¬ 
totic equivalence [H].) Eurthermore, the rapid decay of p allows us to extend the integration range to the full real 
line, so that the integral becomes p{25). Overall, we have the approximation 


H{(jj) ~ 2p{2S)^ (p(cj)e 


iScj\ 


(15) 


provided the errors incurred can be suitably controlled. This is achieved in Appendix]^ with the following result (a 
particular application of Lemma A.l). 

Theorem II.1 Let /3 > 0 and 0 < a < 1. Suppose the function p has transform obeying p{w) = ) and 

p'{uj) = for all w € M, and p''{uj) = for |a;| > ujq, for some ujq > 0. Then the 

integral I{uj) defined in Eq. (12) obeys 


for all Lu, and 


I{uj) = p{2S)p{u;) + 0(|a;|“-ie-'^l‘"l“) 
I'iuj) = p{25)p'{uj) + 0(|a;|“-2e-^l“l“) 


for |w| > ujo- 

Consequently, under the hypotheses of the theorem, we have 

H{uj) = 2p{2S)?R: {p{u})e^^‘^) + 0(|wr"^e-'^l‘^l“) 


(16) 

(17) 

(18) 


for all w, which shows in particular that H{u}) = 0(e The derivative of H will be studied later, and we will also 

give specihc examples of functions p that meet the hypotheses of Theorem m after giving two further refinements 
of our construction. 

Eirst, it will be convenient to work with a function that has a nonnegative transform. This is achieved by the 
expedient of taking the convolution of H with itself, defining a smooth function 


/ OO 

dt' H{t - t')H(t') 

-OO 


(19) 
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which has support [—25,25] and transform K{ui) = H{uj)‘^ obeying 


K{uj) = 2^{25f (|(^(w)|2 + 3? (^(a;)2e2*'5“)) + Odwr-^e-^/^l-r). 


( 20 ) 


The behavior of K{t) near the endpoints t = ±25 of the support is worth studying. Owing to the support of H, the 
integration range in Eq. (19) may be replaced by [—5, <5] and, changing variables to s = t' + (5, we may calculate 


pE 

K{—25±e)= / ds (p{e — s)ip{s)(p(25 — s)(p{25 ± s — e) (p{25)'^ / ds (p{e — s)(p{s) 
Jo Jo 


( 21 ) 


if e ^ (5, using the support of the ip{e — s) factor to restrict the integration range. The last integral may be recognised 
as a half-line convolution and is therefore the inverse Laplace transform of (^(p)^. 

The second refinement arises because K has oscillations of magnitude comparable to its local mean value. It will 
be useful to damp these by a further trick: exploiting again the fact that the derivative of if falls off more rapidly 
than it does, the function 


L{uj) = K{uj) -I- - (^K{uj ± tt/{ 25)) ± k{uj — 7r/(2(5))^ 

has asymptotic form 

L(w) = 4(p(2(5)2|(p(w)|2 + 0(|w|“-ie-2'^l‘"l“) 


( 22 ) 


(23) 


valid for M> ojo for any wq > tt/{25). 
Theorem 11.1 the fact that \uj — 7r/(2(5)|'’'e 


This is proved by using the estimates on ip and p' in the hypotheses of 
-2/3|i^-7r/(2(5)|“ _ ^ following result. 


Lemma II.2 If F{ijj) = /(a;)e^*‘^‘^ where 5 > Q is fixed and f is differentiable, then 


F{u;)±1{F{u;±7t/{25))±F{u,-7t/{25))) = oI sup |/'(C)|). (24) 

^ \K-<^l<7r/(25) / 

Proof Elementary calculation shows that the right-hand side is [/(‘j-') ~ /(^ ~ '^/(2<5)) + /(w) — /(w -I- tt/{25))\ 

and we then use the mean value theorem. □ 

The function L may be recognised as the transform of 

L{t) = {I ± cos{TTt/{25)))K{t) = 2 cos^(7rt/(4^))iL(t), (25) 

which still has support [—25, 25]. Moreover, the cosine factor executes a single period of oscillation on the support of 
K, and so does not disturb the general ‘bump function’ shape. 

We also want to control the derivative L'{u>). Recalling that iL(a;) = 23fi(e*“'*/(a;)), it is easily seen that 

H'{u}) = 253?(fo“'^/(w)) -f 0(|/'(w)|) (26) 


and therefore K — has derivative 


k'{uj) = 2i5 (e*“^/(a;) -b e'^^^r (w)) (w)) + 0(|/(a;)/'(w)|) 

= 4,5tf{ie^^^^I{ujf) ± 0(|/(a;)/'(a;)|). 


Substituting this expression into the derivative of (22) and applying Lemma 11.2 we find 


L'iio) =0[ sup |/(C)/'(C)I = 0(|a;r-ie-2/5for) 


(27) 

(28) 


(29) 


valid for jwj > wq > 7r/(25). 
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B. Specific examples 


A class of specific examples can be obtained as follows. We choose to be the inverse Laplace transform of 
(f{p) = exp(— (pt)“) where r and a are positive constants with 0 < a < 1. The role of r is to set a timescale. As 
usual, the Laplace transform is defined by <f{p) = Jq dt (p(t)e~P*, and the branch cut for fractional powers is taken 
to lie along the negative real axis. Functions of this form have a number of interesting properties, some of which will 
be discussed in Sec. m Given this definition, it is easily seen that (p is supported in the positive half-line: if f < 0 
one may, in the Laplace inversion integral, complete the contour in the right-hand half-plane, in which p is analytic, 
and conclude that ip(t) = 0. Consequently, the Fourier transform tp oi p extends to an analytic function in the lower 
half-plane and may be found as the boundary value 


(p(uj) = lim (p{uj — ie) 
£->■ 0 - 1 - 


lim (pie + ico). 
£-> 0 - 1 - 


This gives 


(p{uj) = exp , 


from which we may read off that 


|(^(a;)| =exp(-/3|a;|“), /3 = r“cos-^. 


(30) 


(31) 


(32) 


The derivatives of p exist except at w = 0, and obey 


W{^)\ 


Q:t“ I i^(cli) I 


W\^)\ < 


/ a{l — a)T°‘ 


l<^H| ■ 


(33) 


In particular, <^"(a;) = 0(|:,c>(a;)| |a;|““^) for uj S [wqjOo) for any wq > 0. The hypotheses of Theorem |IL1| are therefore 
satisfied and the construction above leads to a smooth, even and nonnegative function L : K —>■ [0, oo) with compact 
support in [—2(5, 25] and with Fourier transform that is analytic, even and nonnegative and obeys the estimates 


L(w) = 4(/5(2(5)2e-2/5|‘^r + O 


,-2/3|c^|° 


|l-a 


(34) 


for w ^ 0, and 


-2/3hr 


^ /0 2/3 


( 35 ) 


valid for w > wq with wg > tt/ (2(5). If, as suggested above, 5 is chosen to be the first local maximum of p^ then 5 oe t, 
and (in any case) can clearly be arranged so that the width of the support takes any desired value. 

Meanwhile, the overall integral of L is 1/(0) = -^(0) + K{tt/{25)) = H{0)^ + H{tt/{2S))^, and we may divide by 
this quantity to obtain a function with unit integral all the previously stated properties and additionally having a 
unit integral. Note that the coefficient of in the normalized version of Eq. (34) is then 4p{25)^/{H(0)^ -I- 

H{n/{25)f). ^ 

The approximation (21) gives the behavior in t-space near the end-points, noting that p{p)^ = ip{2^^°‘p) has inverse 
transform 2~^l°‘p(t/2^^^, 


K{-25 + e) 


e ~ 




L{—25 -I- e) ~ -^e^K{—25 + e) 


(36) 


for e <C (5 (here L has not been normalized to unit integral). 

The construction described in this section can be illustrated by taking the case r = 1, a = 1/2, 5 = 1/2 as a 
concrete example. In this case, the Laplace transform <p{p) = e“v^ may be inverted in closed form to give 


pit) 


^-3/2g-l/(4t)y',y4^ t>0 

0 t < 0 


(37) 
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FIG. 1: A plot of H{t) for the case a = 5 — 1/2. 


See, e.g., formula 17.13.31 in Ref. [29]. As S has been taken slightly larger than the first maximum of ip (which occurs 
at (5 = 1/6) the resulting function 


H{t) = ip{t + l/2)p{l/2-t) 


|(l_4i2)-3/2e-l/(l-4t=) |i|<l 

0 |i| > i 


(38) 


has a broad peak with local minimum at t = 0 (Fig. [^. The Fourier transform of H, and its leading order asymptotic 
approximation 


H{uj) 


e ^/^e cos 




(|w| ^ oo) 


(39) 


are plotted in Fig. as the graphs become indistinguishable even at parameters of |a;| ~ 20, they have been plotted 
again in Fig. over a larger parameter range, but with the leading exponential decay factor divided out. The 
function L defined by ( |2^ has integral 0.0658,^ and when one normalizes L so that it has unit integral (and hence 
1/(0) = 1) we obtain the function plotted in Fig. which has L(0) = 1.4990. Near the switch-on and switch-off 
regions, i.e., |t| ^ 1, we have 


L{t) ~ const • 9{1 — |t|)\/l — |t|e 1*1^^ 

The (normalized) transform is plotted in Fig. and obeys 


Liuj) = 2.9324 e 




O 




\/M ) 




(40) 


(41) 


^ The numerical method used was to evaluate 101 samples of L at equally spaced points between 0 and 1 inclusive, using numerical 
integration. The integral of L over [0,1] was evaluated by interpolating a piecewise linear spline through the samples and integrating 
numerically; the result was doubled to obtain the full integral of L. These computations were performed in Maple 18. 









5=0.5 



FIG. 2: A plot of and its leading order approximation for the case a = 5 = 1/2. 


8=0.5 (exponential decay removed) 



FIG. 3: A plot of e“/^ and its leading order approximation for the case a = 5 — 1/2. Here the exponential decay has 

been removed to leave an envelope of approximately constant amplitnde. 
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FIG. 4: A plot of L{t) for the case a = 5 = 1/2. 


5=0.5 



FIG. 5: A plot of L{ui) and its leading order approximation for the case a = 5 = 1/2. Here L has been normalized so that L 
has unit integral. The asymptotic approximation becomes good for |a;| > 10. 
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8=0.5 (exponential decay removed) 



FIG. 6: A plot of e'^ L and its leading order approximation for the case a = <5 = 1/2. The exponential decay has been 
removed. 


which is seen to become quite accurate even at moderate values of |w|. This impression is supplemented by Fig. |^in 
which the same two functions are plotted over a different parameter range and with the leading exponential decay 
divided out - thus the leading order asymptotic contribution in this figure is constant. 

Johnson [3D] has discussed a class of related functions which can be simply defined in coordinate space, and whose 
Fourier transforms obtained by a saddle point approximation have a similar asymptotic form to (39), but with an 
additional power-law decay factor. 


C. An electrical model 


As an example of how a leading edge of the form (40) might arise in controlling a physical device, consider a wire 


with uniform resistance per unit length p and capacitance per unit length c relative to a ground potential. At time t 
and position x along the wire, let I(t,x) be the current flow in the direction of increasing x, and let V(t,x) be the 
potential relative to ground. Elementary application of the laws of Ohm, Kirchhoff and Volta gives 


dV__ dV _ 

l)x~~^’ l)t ~ 

which combine to show that V obeys the diffusion equation 

dV _ 1 d‘^V 
dt pc dx'^ 


1 dl 
c dx' 


(42) 


(43) 


This can be solved for t,x > 0 with specified boundary condition V(t,0+) = Vo(t), representing the applied control 
voltage, and with V(0+,a;) = 0 for all x. Solving by Laplace transform in t and discarding the solution growing in x, 
one hnds 


V{t, x) = f G{t - t'-x)Voit'), G{p- x) = 
^0 


(44) 


Note that G{p;x) = (p{p), with t = x^pc. In a model where the control voltage is applied sharply at t = t*, 
Vo(<) = V/r?(t — t*) and one finds the exact solution 


V(t, x) = 14 erfc 


I x'^ pc 
4(t — 4) 


V.. 


I 4(f 4) 

npcx^ 


(45) 
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for i —>■ where erfc is the complementary error function. Comparison with the approximation (40) shows that 

our specific function L{t) for a = 1/2 gives a plausible model for the rise of a control voltage at some position x along 
the wire. Of course, there are limitations to this model, because one cannot model the switch-off in the same way. 
However, it indicates that the sampling functions used here are less artificial than might be thought at first sight. 


D. Relation to positive stable distributions and Fox H-Functions 


We note some further special properties of the functions used as the basis of our construction. Recall that 
is a function on the positive half-line with Laplace transform exp(—(pr)“). One may show that (/? is a nonnegative 
function with unit integral over IR+ (see, e.g.. Ref. m) — in other words, it is a probability density function. Fixing 
a, the probability density function for the sum of two random variables with densities defined by ti and T 2 is given 
by the half-line convolution of their density functions, which therefore has Laplace transform exp(—(pra)") where 
T3 = (t“ -|- This is another distribution in the same family and for this reason, in conjunction with their 

support on the positive half-line, they are called positive stable distributions. 

We have already given the explicit form of (p{t) in the case a = 1/2; it is worth noting that the distributions for 
other values of a can be expressed in terms of Fox iL-functions. Consider the one parameter family of functions, 
ga(t), whose Laplace transforms are of the form e“^ , 


/ ga{t)e~^^dt = . 

-^0 

A series expansion for gdt) was obtained by Humbert [32] and proved rigorously by Pollard |3T] : 

, , 1 ^ (-1)'= 
gait) = - 

TT 


TT ^—‘ 

k=l 


Tiak 1) sinTrfca . 


Schneider 


has discussed these functions, and shown that 




-1 1 

—a~^ a~^ 


(46) 


(47) 


(48) 


where is a Fox H-function (see, for example. Ref. [Mj). The asymptotic form as t 0'^ is 


where 


g = 


and 


D = 


Ot-^e-^* ", 

(49) 

a 

(50) 

1 — a ’ 

2-a 

(51) 

2(1-0) 

_^)^a/{l-a) ^ 

(52) 

yl/m-a)] 

(53) 

^27r(l — a) 


Thus for the case a = 1/2, we have ly = 1, /j, = 3/2, w = 1/4, and D = l/\/^ so that gi/ 2 (t) has the same switch-on 
behavior as that of (p(t) given in Eq. (37); indeed, gi /2 and ip coincide. For more general values of a, the switch-on 
rate is characterized by the parameters iz and p, which in turn are determined by Eqs. (|50|) and (511. 


At large values of t, it is clear from Eq. (47) that gait) decays slowly - the distribution has a heavy or fat tail, and 
indeed, for 0 < a < 1, pa does not even possess a finite mean. 

Einally, we mention that the functions pa also appear in the consideration of fractional derivative extensions of the 
diffusion equation, formally = Kd^P/dx^. Thus the electrical model of the previous subsection could in 

principle be generalized to a rather broader class of signalling problems in time-fractional diffusion-wave equations 
(see, e.g., [35]) thus modelling the switch-on process by an anomalous diffusion processes. 
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III. HIGH MOMENTS 

In this section, we will give some general formulas for moments of quadratic operators, with special attention to 
the asymptotic forms for high moments. Let T be a normal ordered quadratic operator which has been averaged in 
time with a sampling function. In general, we can expand T as a mode sum 

T = 'y Oi + Bij tti Gj + al Oj), (54) 


where Aij and Bij are symmetric matrices, and a\ and Ui are the usual creation and annihilation operators for mode 
j of a bosonic field. (The expression as a discrete sum is for convenience of presentation.) We define the moments of 
Tby 


Mn = {Tn , (55) 

where the expectation value is in the vacuum state. Thus /tq = 1, Mi =0; formulas may be derived for higher 
moments, such as 


^12 = 2J2B^JB*,, 
-ij 


(56) 


M3 — 4 y ( Bji Aik Bkj , 

ij k 


(57) 


and 


M4 — 3 M2 + y ( (16 Bji Bkt B*j. B*g + 4 Bji Aik Aji Bki + 4 Bji Aik Aki B*^). (58) 

i j k £ 

Consider the example of a massless scalar field, ip. The mode expansion for its time derivative at a fixed spatial 
point can be taken to be 



(59) 


where oj = |k| and C is a quantization volume with periodic boundary conditions. Let T be the time average with 
sampling function f{t) of : so 


In this case, we find 


and 


/ OO 

dtf{t) :<p^ . 

-OO 


^fuJTuTc ^ 

Aje = y f(coj - uje) 




(60) 


(61) 


(62) 


Here, we have defined the Fourier transform as in Eq. and assumed that the sampling function is real and even, 
f{t) = f{—t) = so its Fourier transform is also real and even. 


fioj) = /(-cu) = riuj) . 

We will also assume that / is normalized to have unit integral, so 


(63) 


/( 0 )=/ = 


( 64 ) 
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Because the summands in the expressions for /r„ depend only upon the frequencies of the modes, in the limit of large 
V we have 



In this limit, the factors of V cancel, and the moments become integrals of products of / factors. 

In the case of a Lorentzian sampling function of width r, 

f(t) = fUt) = , J , (66) 

Tr[t^ + T^) 

we have 

/>) = /z,(w)=e-"H. (67) 

If we use this form, the expressions for the /i„ given above reproduce the values given in Table I in Ref. [S]. 


A. The Dominant Contribution to High Moments 


As n increases, the number of terms in expressions for such as Eq. (58) grows rapidly. Fortunately, there seems 


to be one term which gives the dominant contribution. For n > 3, this term is 


M„ = 4 ^ 


R A 4 

^3132 -^3233 -^3334 


■■A, 


-13^ 


B 


3-n3l 


( 68 ) 


while for n = 2 the leading factor is 2. For the case n = 4, this term is the last one in Eq. ( |58[ ). We can numerically 
test the hypothesis that ~ for large n using the case of Lorentzian sampling oi : :, and the data in Table I in 

Ref. [H]. In this case of n = 4, we find M ^/« 0.843. However, this ratio steadily rises toward unity as n increases, 
becoming 0.966 at n = 10 and 0.993 at n = 20. 

We can understand the origin of the above expression as follows: In every term of T” contributing to /r„, there 
must always be a factor of Bj^j^ on the left and one of on the right. This arises because Bij Ui aj is the only 
term in Eq. (54) which does not annihilate the vacuum from the right, and R* aj is the only term which does not 


annihilate it from the left. The overall factor of 4 in Eq. (68) is a combinatorial factor coming from the two ways to 


order the jij 2 indices on the left, and the two ways to order jnji on the right. After the two B factors, there are 
n — 2 remaining factors, which can be combinations of A and B factors. However, the relative minus sign between 


the terms in the argument of / in Eq. (61) tells us that A tends to fall more slowly with increasing uji than does B. 


Thus the largest contribution arises when all of these n — 2 factors are those of A, as in Eq. (68). In addition to the 
contribution M„, which can also be written as a trace = 4Trthere are a number of other terms in 
which the factors in the trace are reordered; these are also suppressed relative to M„. In any case, all of the remaining 
terms in /j,„ are nonnegative, so Eq. (68) gives a lower bound on the moments, M„ < /x,j. 

We may use Eqs. (61) and (62) to write an expression for in terms of Fourier transforms of the sampling 
function. In the limit H > oo, the result will be of the form 


Mn = Cn dhJl--- duJn{uJl ■ ■ ■ W„)^ /(wi + UJ2) f {uJ2 - UJ3) ■ ■ ■ /(Wn -1 “ UJn)f{^n + Wl) . 


(69) 


Here 


Cr,. = 


(27r^ 


(70) 


and p = 3 for the case of given in Eq. (60), which will be our primary example. Equation (69) can also be obtained 
as one of the terms appearing in the evaluation of p,„ by Wick’s theorem, as in Sec. HI.A and Appendix B of [5], 
where it was identified as the dominant contribution in the case of Lorentzian smearing. 
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B. Asymptotic Form of the Moments 


In this subsection, we develop a procedure for finding asymptotic limits of M„ for large n. Consider the case where 
n is even, so we may write n = 2m and rewrite the integral in Eq. (69) as 

/*C30 

M2m = C2m / duJiduJm+1 (oJlLOm+l)^ {Gm-l (wijWm+l)]^ , (71) 

JO 

where Gm-i is defined by an m — 1 dimensional integral, 

poo 

G'm-l(wi, Wm+l) = / dljJ 2 • ■ ■ du!m { 0 J 2 ■ ■ ■ OJm)^ f {OJI + UJ 2 )f{ 0 J 2 — OJ 3 ) • ■ ■/{oJm — OJm+l) ■ (72) 

Jo 

Note that in general, Gm-i will not be symmetric in its two arguments. We can also write an expression for a general 
value of n as 


M 


m+m' — 


poo 

/ dwdfl Gm-l(w, n) Gm'-l(w, fl) . 

Jo 


The definition, Eq. (72), leads to a recurrence relation 

poo 

G,iu;,n)= dU^f(^-OGi-i{uj,0, 


(73) 


(74) 


with Go(w, fl) = f(uj + ^). 

Thus, knowledge of Gi allows us to compute all of the G^, and hence the dominant contribution to the moments, 
M„, for all n. 

We are especially interested in the case of compactly supported sampling functions, such as those discussed in 
Sect. 1^ If we adopt units in which r = I, the Fourier transform, /(w), will fall faster than any power of to when 
w ^ 1. We will assume that it falls as an exponential of a;“, where 0 < a < 1, specifically that 


/(w) = ye 




o 




ll-a 


(75) 


for w 7 ^ 0. We are interested in the asymptotic forms of Gi, Gm, and finally of M„. To start, let us proceed 
heuristically. As / becomes more flat at large values of its argument, a reasonable conjecture is that 

Gi(w,fI)~27r/(0)fIP/(a; + fI). (76) 

Here /(O), the sampling function evaluated at t = 0, is expected to be of order unity. In making this approximation, 
it is important that our the assumption (75) entails that / does not oscillate at leading order (cf. the behaviour of 
the function H in Eq. (dU)); if it did, one might expect a similar approximation, but modified by a factor which 
would feed through to later estimates. Proceeding with the conjecture (76), we may now use the recurrence relation, 
Eq. (74), to find a corresponding conjecture for the asymptotic form for Gm- We make the ansatz that 

G^) (cu, n) cs [2^/(0)]"* f(uj + H). (77) 

(78) 

(79) 


where we have added a superscript to denote the value of p. Insertion of this ansatz in Eq. (74) leads to 

pOO 

Gtl^ [27rf{0)r / dC /(II - /) /(w + /) = [2^/(0)]™ (cc, H). 


Now the asymptotic form, Eq. (76), leads to 

'lip) 


G^Zh (w. ^ {27rf(0)r+^ /(cc + H). 


which proves our ansatz by induction at this (currently heuristic) level. 

In Appendix [B| these conjectures will be established rigorously. To be precise, it will be shown that for each p > 1 

'—' (n) (p) 

and TO > 1 there exist polynomials Pm and Qm of degrees mp — 1 and mp respectively, with nonnegative coefficients 
independent of uj and H, such that 




ip) 


(n) 


(a; + 


/3(u;+Q)° 


(80) 
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for all WjH > 0. This shows that the error terms in the approximation are suppressed relative to the dominant 
term. For instance, if —>■ oo, with w/fl fixed, the error term is suppressed by a factor 11““^ relative to the main term. 
The statement (80) also holds if m = 0, under the convention that a polynomial of degree —1 vanishes identically; 
in this case we simply reproduce the assumption (75) on /. Although we do not control the polynomials Pm ^ and 
Qm'^ explicitly (for example, to show that the coefficients are independent of m), we do show that there are positive 
constants A, B, c and D such that 


^^m^mpe-fo+o)“ < G!^\uj,n) < {cD"^n^P + 


(81) 


for all w, n > 0 , where Rm^ is a polynomial of degree mp — 1. These rough bounds provide further evidence that ([77]) 
is robust. 

Now we seek the asymptotic form for M„ when n ^ 1. Use the asymptotic form for Gm, Eq. ( [77| , in Eq. ( [73| to 
write 


M„ ~ Cn [27r/(o)]"-2 / dujdnujpn^^-'^'fp p{uj + n), 

Jo 

where n = m + m'. Next let u = a; + U, and write 

1*00 pu 

M„ ~ [27r/(o)]"-2 / dup{u) / dn {u - n)p n^^-^'>p, 

Jo Jo 

and then let x = Q/u to obtain 

pOO pi 

M„ ~ [ 27 r/( 0 )]"-^ / dup{u)u'^P+^ / dx{l - x)p . 

Jo Jo 

The integration on x may be performed to yield 


~ Cn [27r/(0)] 


^_2 p'-[in-l)p\l r 

{np+iy. Jo 


dup{u) u^P+^. 


(82) 


(83) 


(84) 


(85) 


If n 3> 1, the factor of will cause the dominant contribution to the u-integration to come from the region where 

M ^ 1, so we may use the asymptotic form, Eq. ([7^ to find 


~ Cn [27r/(0)] 


j -2 p'\{n-l)p]\ 

{np+iy. a (2,5)(”P+2)/° 


{np + 2 ) 
a 


( 86 ) 


The most important part of this expression is the final gamma function factor, which leads a rapid rate of growth of 
the high moments, M„ oc (pn/a)!. This result will be used in Sec. IV to infer the tail of the probability distribution. 
Although we have not performed rigorous error estimates here, they could be done; however, the rough bounds (81) 
show that the gamma function factor will be a robust feature of the result. 

The argument given here and in Appendix does not apply to the case of a Lorentzian sampling function, where 
a = 1, as the error term in Eq. will now be as large as the dominant contribution. A different argument specifically 
for the Lorentzian case is given in Appendix This argument reproduces the asymptotic form of the moments found 
in Ref. [5]. 


C. Numerical Tests 


Here we discuss some numerical tests of the asymptotic forms derived in the previous subsection, focussing in 
particular on the form of Gi{uj, H) for large H with fixed w/H, given in Eq. (76). Here we wish to test this result in 


the case that /(w) is exactly given by Eq. (751, with /? = 7 = 1, so 

/(w) = e-l-l“ . 

Now the exact form for Gi{u!, ft) may be expressed as 

Gi{uj,U) = J dz/z/P exp 


(87) 


l-z/l 




( 88 ) 
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R 



FIG. 7: Here R, the ratio of the numerical integration for Gi, from Eq. 
plotted for the case a = 1/2 as a function of fl for various values of w/fl. 


to the asymptotic form, Eq. (891, for p = 3, is 


R 



FIG. 8: A plot of R for the case a = 1/4. 


where = ^/n. If we use/(O) = (27r) ^dw/(w) = r(l/a)/(7ra), then the asymptotic form in Eq. (76)maybe 
written as 


OP 

GiA{uj,n) = —r{l/a) e-(‘"+^)“. 
a 


(89) 


Let R be the ratio of the result of numerical integration of Eq. (88) to the asymptotic form, Eq. (89), for p = 3. 
In Fig. this ratio is plotted for the case a = 1/2 as a function of fl for various values of uj/fl. We see that i? —)■ 1 
for n > 10^, as expected. (For the case w/fl = 0.1, this has not quite occurred on the scale plotted, but does occur 
for larger values of ft.) This plot is repeated in Fig. for the case a = 1/4, where the asymptotic limit is reached for 
^ ^ 10^- The case a = 1/8 is illustrated in Fig. ^ where the asymptotic form holds for ft > 10^°. Note that before 
the asymptotic limit is reached, Gi(w,r2) tends to exceed the asymptotic form, in some cases by large factors. This 
emphasizes that the asymptotic estimates are lower bounds on the exact results. Note that i? —>■ 1 for large 12 in all of 
the cases plotted, including when a < 1/4. The only case plotted where i? < 1 in part of the range, is in Fig.[^ when 
a; < n. However, this effect seems to be small, and does not alter the conclusion that the asymptotic form, Eq. ( |8^ , 
holds for large H. 

A separate numerical test can be obtained by numerical evaluation of from Eq. (69) using the form of / in 


Eq. (87). Unfortunately, this is only feasible for relatively small values of n. For the case of p = 3 and a = 1/2, this 
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R 



FIG. 9: A plot of R for the case a = 1/8. 


evaluation for leads to M 4 « 3.0 x 10^^ 
to M 4 « 3.0 X 10^2 and Mq « 1.0 X 10^^ 
small to expect detailed agreement. 


and Mq « 1.5 x 10^®. By comparison, the asymptotic form, Eq. ( 861 , leads 
These results are in fair agreement, despite the fact that n = 6 may be too 


IV. THE TAIL OF THE PROBABILITY DISTRIBUTION 


In this section, we will use the asymptotic form of the moments, Eq. ( 861 , to obtain information about the probability 
for large fluctuations. Let the probability distribution, P(x ), be a function of the dimensionless variable x = T 
Our main example is p = 3, so a; = Tr^ 

using units in which r = 1. We can now restore an arbitrary value for r by using the variable x. 


we computed the asymptotic form for the moments 

We follow the 


discussion in Sec. V of Ref. [B], and postulate that the asymptotic form for P{x) for a: 1 can be written as 


P{x) 


Co x'^ 


(90) 


for some constants cq, a, b, and c. The n-moment for n ^ 1 is given by 


fJ-n = 


/ oo pc 

x^ P{x) dx ~ Co / 

-Xn J 0 


^n+b 


_ G) ^_(n+b+l)/c 


r[(n + b+ l)/c]. 


(91) 


We assume that the asymptotic form of the dominant contribution to /i„, given by Eq. ( 86 ) will match Eq. (piT 


allowing us to infer values for the constants in Eq. (90). Note that the asymptotic form of the moments is growing 


too rapidly with increasing n for the moments to define a unique probability distribution (at least according to 
the Hamburger or Stieltjes criteria m)- However, the difference between two probability distributions with the same 
moments must be an oscillatory function, all of whose moments vanish. We are less interested in the detailed functional 
form of P{x) than in the its average rate of decay for large x, and we expect this average rate to be of the form of 
Eq. (90). This reasoning can be supported by rigorous arguments - see of Sec. VI of Ref. [ 6 ]. 

First, note that 


np - 


and that 


= {np/a + — I)!, 


v-(p+i) 


(92) 


[(n- l)p]! ^ _ 

{np + I)! {np/a + l/a){np/a) ■ ■ ■ {np/a — 2^) 

Q,-(p+i) 


{np/a + 2/a — I) • • • {np/a + 2/a — p — 1) 


1 + 0 I - 

n 


( 93 ) 
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In the last step, we use the fact that n is large, but p is fixed, so 


1 


1 


1 + 0 ( - 
n 


np + A np 

for any constant A. We assume that Cn is of the form 

0 „ = Bo 


for constants Bq and B. This form holds for the case of 93 ^, given in Eq. (70|. 


Next we drop the 0(l/n) terms in Eq. (93) and rewrite M„ in Eq. ( [ 86 | ) as 

M„ ~ Bop\ 7 ^ B^ [27r/(0)]”-2 ( 2 / 3 )-("p+ 2 )/“ r{np/a + 2/a-p-l). 


Now we can compare Eqs. (911 and (96) to obtain expressions for b and c in terms of p and a alone: 


a 

c = — 
P 


(94) 


(95) 


(96) 


(97) 


and 


b = c ( ——p — — 1 = -— ^ — (a + 1 ) . 

a J p 


We also obtain expressions for a and cq as 


a = 2/3 [2^/(0)i3]-“/P, 


(98) 


(99) 


and 


Co = ca^’^+^^/^Boph^ (2/3)-(2/“) [27r/(0)]- 


( 100 ) 


We are especially interested in operators, such as stress tensors, for which p = 3, so c = a/3 and b = —{Aa + l)/3. 
Although Eq. ( | 86 [ ) does not strictly hold for the Lorentizian sampling function, a = 1, we do reproduce the c = 1/3 
result of Ref. [ 6 p Compactly supported sampling functions with a < 1 lead to an even slower rate of decrease of the 
tail, such as c =1/6 for a = 1/2. Recall that the asymptotic form for the high moments which we have used, Eq. ( [ 86 | , 
is expected to be a lower bound on the exact moments. This means that the form of the asymptotic probability 
distribution, Eq. (90) is a lower bound on the exact P{x), which could potentially fall even more slowly for large x. 


V. SUMMARY AND DISCUSSION 


In this paper, we have studied the probability distributions for quantum stress tensor operators averaged in time 
with functions of compact support, which describe measurements in a finite time interval. Our main interest is in 
the asymptotic form of the distribution, which give the probability for large fluctuations in the vacuum state. This 
asymptotic form is determined by the rate of growth of the high moments of the distribution. This rate of growth 
was inferred by first identifying a dominant contribution to the high moments, Eq. ( 68 ), which is determined by 


the Fourier transform of the sampling function. We next used asymptotic analysis to find the the rate of growth of 
the dominant contribution, which is given in Eq. ( |96[ ). This rate is typically too fast for the moments to uniquely 
determine the probability distribution according to the Hamburger or Stieltjes criteria. However, asymptotic form for 
the moments is adequate to give estimates for the probability for large fluctuations, which was done in Sec. |IV] Note 
that the dominant contribution, M„, is a lower bound on the exact moments, /r„, so our estimate of the asymptotic 
probability distribution is a lower bound on the exact result and hence on the probability of large fluctuations. 

Compactly supported sampling functions have Fourier transforms which fall more slowly than an exponential, 
resulting in more rapid growth of the moments compared to the case of Lorentzian averaging, and a more slowly 
falling tail of the probability distribution. For example, the sampling functions with a = 1/2, such as those discussed 
, result in moments growing as ( 6 n)! and an asymptotic probability distribution for large x of the form of 
where x is the dimensionless measure of the sampled stress tensor. This slow rate of decrease leads to the 
possibility of significant physical effects from large fluctuations. Furthermore, the rate of decrease is very sensitive to 
the value of a, with smaller a leading to an even greater probability of large fluctuations. The value of a is linked to 
the rate of switch-on and switch-off of the sampling functions, as described by Eqs. (49), (50), and (51). 


in Sec. H B 
^-1 


w 
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There are several possible physical applications of these results, including the effects of quantum radiation pressure 
fluctuations on barrier penetration by charged particles or upon lightcone fluctuation effects in a nonlinear material. 
These are quadratic operator fluctuation effects, analogous to the linear electric field fluctuations studied in Refs. mm- 
In both cases, the physical system may define a choice of the sampling function /r(t). For example, the effects of the 
fluctuations of the squared electric field in a nonlinear material are determined by an integral along the path of a probe 
pulse, as described in Ref. [3]. This integral is in turn depends upon the density profile of the slab. This illustrates 
the general principle that physical processes which begin and end at finite times should be described by sampling 
functions with compact support, and the specific form of these functions is ultimately determined by the details 
of the physical process in question. Because all of the probability distributions treated here fall more slowly than 
exponentially, even at finite temperature sufficiently large fluctuations of, for example, energy density are more likely 
to be vacuum fluctuations than thermal fluctuations. This and other physical consequences of the results described 
in this paper need to be more carefully explored. 
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Appendix A: Integral estimates 


The aim of this appendix is to prove Theorem |II.1[ This is accomplished as a special case of the following result. 

Lemma A.l Suppose that positive functions 4','!']^ and T are defined on (0, oo), with iki monotone decreasing, and 
obeying 

(Ja;(l + a;)T(a;) < oo (Al) 

and 

poo pXuj 

4'(w) / doj'T{uj')=0{^i{oj)), / da;'w'4'i(w-a;')T(w') = 0(4-1 (w)) (A2) 

J\iij J 0 

foruj > 0 andfixedX G [0,1]. Lett/) : (0,oo) —>■ C be differentiable andv : K —>■ C fee continuous, withv{uj) = 0(T(|a;|)) 
for all uj gM. and suppose there is wq > 0 so that tpitv) = 0(4'(a;)), if'iuj) = 0(4'i(w)) for uj > ujq. Then 




duj' tfiio + ui')v{ui') = if>(uS) 



dio' v{u}') + 0(4'i(a;)) 


(A3) 


for u! > uio/{l — X) if X < 1 or UJ > 0 if ujQ — 0. If, additionally, if' is dif feren tiable with if"{uj) = 0(4-2 (w)) where 4-2 
is a positive, monotone decreasing function on (0,oo) so that properties (A2) also hold with 4- and 4-i replaced by 4-i 
and 4-2 respectively, then 


d 

dui 



duj' if{uj + uj')v{uj') = if'{uj) 



duj' v{uj') — A'0((1 — A)u;)r;(—Aw) + 0 ( 4 - 2 (w)) 


(A4) 


holds for u! > uJo/{l — X) «/ A < 1 or w > 0 z/wq = 0. 


Proof: The integrals in Eq. (A3) are finite because v is necessarily absolutely integrable and if is bounded. A 


straightforward rearrangement of terms gives 

/ OO poo p — XuJ poo 

duj'if{uj + uj')v{uj') — if{u!) / duj'v{uj') =—if{uj) / duj'v{uj')+ / duj'{if{uj + uj') — if{uj))v{uj') (AS) 
-Xlj J —oo J —oo J —Xuj 


and the first term on the right-hand side is 0(4-i(|w|)) by the hypotheses. Next, the mean value theorem gives 
if{uj J- uj') — if{uj) = uj'if'{f) for some C lying between ui and w -I- uj' , giving 

\if{uj J- uj') — ifiuj)\ < Iw'l 4-i(min{a;,a; -I- uj'}) 


(A 6 ) 
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on recalling that 'I'l bounds "0' and is monotone decreasing. Splitting the integration range at w = 0, this gives 


dui' (■i/'(w + — tp{oj)) v(oj') 


' —Xu) 


nOO pXoJ 

<'4'i(a;) / duj'uj'T{uj') + / da;'a;'d'i(a; — a;')T(w') (A-7) 

Jo Jo 


where we have also changed oj' —>■ —w' in the last integral. As both terms on the right-hand side are 0(4'i(|a;|)) the 
first part of the result is proved. For the second part, we first note that 


/ OO nOO 

dui'tf}{u} + uj')v{uj') =— X)u!)v{—Xu!) + / duj'tjj'{oj + ijj')v{w ) 

-Xuj J —Xuj 


d 

duj 


(AS) 


and apply the first part of the lemma to the second term. □ 

The proof of Theorem|II.l|now amounts to verifying that the hypotheses of Lemma [A. 1| are satisfied for the functions 


d'(w) = T(tLi) = exp(—/3a;“), d>i(tLi) = 


exp(—/3a;“) 


,1-a 


«' 2 (w) = 


exp(—/3w“) 


,2—a 


where /3 > 0, 0 < a < 1, in the case A = 1/2. Condition (Al) is obvious, while we may also compute 

^l-ag-/3n“ 




a/3 


{n 




(A9) 


(AlO) 


using Laplace’s method (see Lemma [A. 2| below). As the right-hand side is bounded for all fl > 0 and decays rapidly 
as n — OO, it follows that, for each 7 > 0, the function oj'^ dcj' exp(—/3(uj')°‘) is bounded for all w > 0 (with a 
bound depen ding on 7 ). Hence the first condition in ( |A2[ ) holds for T and d>i and also for d'l and '^' 2 - The second 
condition in (A 2 ) is verified by noting that 


fUj/2 


dui 


{(jJ — U}')'^ a 


r-2 “ ,,2/a-l ,, 


(1 - uV“)’; 


(All) 


where we have ma de th e change of variables v = {uj'/uj)°‘. The integral may be estimated by Laplace’s method in 
the form of Lemma A.2 below; in particular, the integral is 0(w“’^e“^“°) for w > 0. This calculation shows that the 
second condition in (A2) is verified for both d'l and 'I '2 for 0 < a < 1. Note, however, that this condition would fail 
in the case a = 1. Eq. (16) in Theorem 11.1 is now immediate, while Eq. 0 follows on noting that 


d'(a;/ 2 )^ = e 


2 _ „-2/3((^/2)“ _ 


= 0(w 




) 


(A12) 


as 0 < a < 1 . 

It remains to state Laplace’s method in the form we use it. 


Lemma A.2 Let y ,7 > 0 be constant and suppose g : [0, E] —)■ K is continuous with g{0) ^ 0 and h : [0, E] —>■ M 
is continuous, attains its maximum only at v = 0, and has a strictly negative one-sided derivative h'{Q). Then the 
integral 


obeys 


m = f 

Jo 


dv 


,7-1 


ff( 


v}e 


th{v) 


(A13) 


r(7)g(0)e*^(°) 


(t —> 00 ) 


(A14) 


and I{t) = 0{t T'e*^(°)) on (0, 00 ). The asymptotic equivalence (A14) also holds [0, E] is replaced by [0, 00 ), under 


the additional assumptions that h(v) 
I{t) = 0 (t“'>'e*^^°)) on {to, 00 ). 


as V ^ OO and that I{t) converges for t = to; 


this 


have 


Proof The asymptotic equivalence holds by a straightforward adaptation of the argument in §4.3 of [28] (which deals 
with the 7=1 case). As I{t) is evidently bounded on any compact set, we may deduce that I{t) = 0{t~'^e*^^°^) on 
( 0 , 00 ) or {to, 00 ) as appropriate. □ 
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Appendix B: Asymptotic behavior of G-n 


The aim of this appendix is to prove statement (80) on the behavior of Throughout, we will assume 

(75), that /(w) = + O where 0 < a < 1 and / 3,7 > 0. It will also be useful to note that 

the magnitude of f{^) satisfies 


|/(a;)|<ce-^l‘"l“. 


(Bl) 


for some constant c. This follows because / is bounded and obeys Eq. (75), although in general c 7 . 


1. Preliminary estimates 


Define a family of integrals 


i:{oj,n)= / de 


(cj + 




(B2) 


where w,D>0, 0<i^<l and r > 0 is an integer. Noting that the factor (w + fl + ^) ’^e /3(“+f^+C)“ is decreasing on 
[ 0 ,oo), we obtain the obvious estimate 


j:(cc,d)< 


/ 3 (ca;+S 1 )“ 

(lo + Q)’' 




where we dehne, for q > 0 , 


POO 

Jo 


(B3) 


(B4) 


which is evidently a polynomial of degree r in D with nonnegative coefficients (independent of w) that can be evaluated 
in terms of Euler E-functions: for instance, one easily sees that 




r((g+ l)/a) 


D’’ -|- lower order terms. 


(B5) 


In particular, the leading coefficient is independent of r. 

Next, let us make the elementary observation that the inequality 


x°' + y°‘ > {x + j/)“ -f (1 — a) min{x, y}°‘ (B 6 ) 

holds for all a:, ?/ > 0, 0 < a < I. This is seen by the calculation (assuming without loss of generality that x <y) 

x°‘+y°‘-{x + yY = x°‘-a f dz 2 “-^ > x°‘ - axy^-^ > (I - a)x“ , (B7) 

dy 

where the integral was estimated by mutiplying the measure of the integration range by the maximum value of the 
integrand; note that x°‘ < y°‘ but y°‘~^ < 

Now consider 


ft) — f df 

" ’~Jo ^(cc + D-O^' 

where a;,n>0, 0<z^<I, q>0 and r > is an integer. Evidently 


Ky{u;,n)< 


D’’ 


< 


(w -I- Jq 

^r^-P(uj+n)° i-n 


{ui + fl)'^ Jq 




(B 8 ) 


(B9) 

(BIO) 











22 


where we have used the fact that {fl — /{uj + 12 — Y)'' is decreasing on [0,12] because Q < v < r. The last integral 
can be rewritten 


^min{n, ^ (a;+Q)} 






’ min{r2, ^ (cj+f2)} 




9p-d(l-a)(0+w-{)° 


(Bll) 


the first term of which can be bounded by a constant independent of oj 
term, we may estimate the exponential from above by 


replacing the upper limit by oo). In the second 
whereupon the integral is bounded above by 


12'^+^e “)(^/2) and hence by a constant independent of w. Thus the integral in (BIOI is bounded by a constant 


and we have shown that 


Ky{^.^) = 0 


f|rg-/3(c^+n)“ \ 

(w + 12 )'^ J 


(B12) 


for w,12 > 0. Note that the controlling constant in this estimate (i.e., the supremum over 12 > 0 of the integral in 
formula ( B10| )) is independent of r. 

A further family that will appear below is defined by 


Li{uj,n)= de 


(w + O"' 


,-/3(|n-€r+|<^+4r) 


(B13) 


where, as before, a;,12>0, 0<^<1 and r > is an integer. By changing variable ^ i—^ — 12 and then splitting the 
integration region into [0, oo) and [—12, 0], we find (changing variables again ^ i—)■ —^ in the second of these) 


Li{uj,n) = rYuj,n) + KYY^,^) 


and hence deduce that 


Li{uj,n)< 


P(12)e 


—/3(a;+r2)° 


{uj + fiy 


(B14) 


(B15) 


where P is a polynomial of degree r with nonnegative coefficients independent of uj and 12, and whose leading order 
coefficient is independent of r. 


2. Rough bounds on Gm 

The iterative scheme for computing the functions Gm may be written in the form 

= /(^ + . (B16) 

where we have inserted a superscript to track the value of p > 1 and the integral operator is defined by 

POO 

(s(p)G)(a;,12)= / - 0Giuj,0 ■ (B17) 

-'0 


Note that, due to the assumption (B11, we have an estimate 


poo 

|(s(p)G)(cc,12 )| < c / ■ 

Jo 


As a preparation for the main estimate, we first show that 

|GW(c.,12)|<QW(12)e-^(“+^)“ 


(B18) 


(B19) 


(p) 

for w,12 > 0, where Qm is a polynomial of degree mp with nonnegative coefficients independent of ui and 12. The 
statement is true by de finition for m = 0, so suppose it holds for some m > 0. Using the inductive hypothesis and 
also the estimate (B18) gives 


poo 

|g|^^i(cc,12 )| < c / 

Jo 


(B20) 
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which can be expanded as a linear combination of integrals of the form Lg(uj,i2) with p < r < (m + l)p, in which 
the coefficients are positive and independent of a;,n. Using (B15), we deduce that r2)| < 

where is a polynomial of degree (m + l)p with nonnegative coefficients independent of w and fl. Thus the 

bound (B19| holds for all integers m > 0 by induction. 

In principle, the polynomials Qm could be computed explicitly, but this would be a matter of some tedium. The 
leading order coefficient can seen to grow no faster than geometrically with to, because the leading order contribution 
to the right-hand side of inequality (B20) is a constant multiple of which is in turn estimated by an expression 
in which the highest order power of has a coefficient independent of to and p. Accordingly, we have a bound of the 
form 


(w, n) < -h lower order) . (B21) 

To complement the rough upper bound ( |B19 ), we now develop a lower bound. For this purpose, we suppose that 

/(w) > Ae-^l"'l“ (B22) 

for some A > 0. We will show that 

G^) (w, n) > . (B23) 


where B = Ar(l -|- a“^)/(2/3)^/“. By our assumption (B22), this statement holds for to = 0, so suppose it is also 
true for some to > 0. Using the definition of and the fact that all the functions appearing in the integration are 
positive. 



pOO 

/ g-/3(if2-cr+i^+cr) 

Jo 

(B24) 

> A^B'^ 

poo 

/ d^iQ + ^)(Bi+i)p^-f}{G+G+nHr) 

Jo 

(B25) 

> A^B'^ 

pOO 

^(m+l)p / ^^g-/3(^“+(GJ+n+^)“) 

(B26) 


where we have first discarded the integration over [0, fl] and made a change of variables ^ i—> ^-1-12, and then discarded 
all but the highest power of fl. Now 


/*oo /*oo 

Jo Jo 


r(i + »-‘) e,. ^ B e,. 

(2,5)'/" A 


for all cr > 0, using {x + y)°‘ < x°‘ + y°‘. Thus 




(B27) 


(B28) 


establishing the inductive step. 


3. Refined estimate of Gm 


We now wish to prove our statement (801 on 
any to > 0, let TZm'^ be the set of functions H : 


the behavior of G^m , for which we require some further notation. For 
(0, oo) X (0, oo) —>■ C obeying a bound of the form 




Pin)- 


Qin) 


(w -I- 


^—P(uj+Q)° 


(B29) 


for all w, n > 0, where P and Q are polynomials of degrees mp — 1 and mp respectively, with nonnegative coefficients 
independent of uj and n (but depending on the function H). Here we apply the convention that a polynomial of 
degree —1 is identically zero. By the triangle inequality, TZm^ is closed under taking sums of functions, and it is also 

(p) 

closed under constant scalar multiples, so TZm forms a vector space in this way. We also note, using the results of 
Appendix |B H that 


/g, , Lq G for 0 < r < mp — 1, g > 0 


(B30) 
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and 


Ll_^ G ior 1 <r <mp, q> 0. 


It is also clear that if Hi G TZm^ and H 2 {uj, il) < i?i(w, O), then H 2 G TZ 
Given these dehnitions, (80) can be reformulated as the assertion 


(p) 

m • 


for all TO > 0 and p > 1, where 


n(p) _ pip) a 7?(p) 

^ m ^ 

F^Hu;, n) = [27r/(0)]™ /(cc + fl). 


It will be convenient to prove an equivalent version of (B32), namely that 


nip) _ nip) fz T? ip) 

^md ^ '^m 


for all TO > 0, where 


= [ 2 ^/( 0 )]"* . 


(B31) 


(B32) 


(B33) 


(B34) 


(B35) 


The equivalence of statements (|B32|) and (|B34|) is due to the fact that fif ^ G Ti[$ as a result of Eq. (l7^ , and 

iv) 

because TZm is closed under taking sums of functions. _ _ 


We now give an inductive proof of assertion (B34). As (B32) holds trivially for to = 0, the same is true of (B34|. 
So now suppose that (B34) holds for some to > 0. Recalling that = E^P'>Gm\ we may write 


Mp) _ nip) 

'-'m+l '-'m+li 


t(p) f nip) _ nip) 'l _i_ /"7 t(p)^(p) _ nip) \ 

y-'m '-'mdj ^ m+1 d) ' 


(B36) 


The inductive step is established if we show that both terms on the right-hand side belong to TZ^. 

a. First term By the inductive hypothesis, G^ — G^^^ G 7Zm\ so our task is to prove that E^P'>TZ^m C 
i.e., maps functions in TZm^ to functions in TZ^\_i. To do this, it is enough to show that maps every function 
n) = n’’(a; -I- g-^(w-i-n) element of TZ^]^_i, in the cases 0 < r < mp — 1 and = 0, or 0 < r < mp 

and h' = 1 — a. Using the bound (B18), 


(B37) 


poo pp-\-r 

\{E^P^H:){ 0 J,n)\ < c / \ e-^(l^-gl°+l^+-l°) =cLg+’~(a;,0) 

Jo Is + 

and therefore our observations (B30) and ( |B31[ ) prove that E^p'^HI G TZ^m+i required ranges of r and v. 

b. Second term We start by rewriting through a change of variable as 

/ OO 

d^ (e + r!)(-+i)P f{-0 . (B38) 

-n 

Next, using the fact that f{~0 = 27r/(0), we write the above expression as the sum of three terms 

(sWG^),)(a;, n) = G^Zii di^, , (B39) 

where 


Ri (ca, U) = -G^^),(cc, Q) Qp I d^fi-^ , 


n-O 


(B40) 


and 


R2iu:, n) = [ 2 ^/( 0 )]™ 7 / fi -0 + 0 - dim 

J-n 


(B41) 
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where 




(B42) 


Here 

'-'m+ld 
)(p) 


(a;,f2) is expected to be the dominant term, and we need to show that the remainders, i?i and i? 2 , lie 
in and are therefore subdominant. The magnitude of Ri is determined from Eq. (AlO), which reveals that Ri 

is suppressed compared to Gmd(w,H) by a factor of order as H —)• oo, and of order OP for all H > 0. 

Combining these estimates, we may bound |i?i(a;, H)| by a constant multiple of G^^^(a;, H) and conclude that it is in 


Tt. 


(p) 

771+1 ■ 


To control i? 2 , we apply the mean value theorem to g, obtaining 

|5(H + 0 - gm = + ,7)1 < Id ({m+ l)p(H + r 7 )(-+i)p-i + 


a/3(H + ,7)(™+i)p 

(w + H + 77 )!““ 


(B43) 


where 77 lies between 0 and d (If one differentiates g, the two terms in parentheses appear with a relative sign, but 
we have changed this to an addition for the purposes of obtaining an upper bound.) If ^ < 0, we have ^ < 77 < 0 and 
therefore 

|5(II + 6 - 5(11)1 < Id (^(m + , (B44) 

using the fact that (H + 77)l'"+^l?'/(a; +11 + 77)^““ is increasing in 77 on [d 0 ] because (771+ l)p > 1 — a > 0 . Accordingly, 
the contribution to R 2 from [—H, 0 ] may be estimated by 

( n//90("*+l)P \ 

(m + n). (B45) 

On the other hand, if ^ > 0, we have 0 < 77 < ^ and therefore 

|5(0 + d - 5(11)1 < (^(m + l)p(0 + , (B46) 

so the contribution to R 2 from [ 0 ,oo) can be estimated as 

|i?26| < [2^/(0)]"^ 70 ((m + 1)P (b47) 


iv) 

We see immediate ly th at both i? 2 a and i? 2 h belong to 
hence establishes (B32| for all tti > 0. 


which completes the proof of the inductive step and 


Appendix C: Asymptotics for the moments of a Lorentzian 


In this appendix, we perform an asymptotic analysis of the growth rate of the moments for the special case of a 
Lorentzian sampling function. In this case, the Fourier transform of the sampling function, in r = I units, is given 
by Eq. (67). Now Eq. (74) may be written in the case f = 1 as 


Gi{uj,fl) = J exp fl ^^+77 ^+|I —i/| I 


(Cl) 


We define Gia(w, H) as the contribution from the interval 1 < v < 00 , and change the integration variable to s = 1, 

leading to 

r°° ^ 

Giaicu, n) = HP+i / ds (1 + s)P e-^(2s+i+c^/n) ~ op g-G+o) ^ (<^2) 

Jo 2 


In the last step, we used (1 + s)p « 1, as the region where s 1 dominates for large H. If we define Gib(a;,H) as 
the contribution to Eq. ( |C1| ) from the interval 0 < < 1, we may use oj/^l + u + \1 — v\ = a;/H + l to write an exact 
result: 


Gib{oj, H) 


p + 1 


g—(cj+n) 


(C3) 
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Now we have the asymptotic form 


Giiuj, ri) — Gia(a;, $2) + Gihiyi, f2) 


p + 1 






This result can be compared with one obtained from Eq. (A4) of Ref. [6], which leads to 




GiSi) = ^ E 


r=0 


(2f2)’' 

r! 


for which the r = p + 1 and r = p terms reproduce Eq. (C4) of the present paper. 
The recurrence relation, Eq. (74), may be expressed for the Lorentzian case as 

pOO 

^0 


In this case, we replace Eq. (0 by the ansatz 

G^P^UJ, n) « 

where 




«! = 


P + 1 


and bi = 1. 


Next we use Eq. (C7) in Eq. (C6) to write 


G^Zhi^. = am ^) + Y f2) 


^ j7(’^+i)(p+i)e“(‘^+^) 
which leads to the recurrence relations 


. 1 

_(to + 1)(p+ 1) 2r2 


^m+1 — 




(to + l)(p + 1) - 1/J 


(to + 1)(p + 1) ’ 


and 


^m+l — 4” 


(to + l)(p + 1) - 1 ’ 


These relations may be combined with Eq. (C8) to find 


and 


bm, — 


“ to! (p + 1)"* ’ 

(to + l)(p + 1) - 1 
(to — 1)1 (p + l)™“i (2p + 1) 


We may now combine Eqs. (73) and (C7) to write an asymptotic form for the even moments as 


POO 

^2m = C2m / dujdO. \Grn-l{^ 
^0 


= G, 


p! [(2 to - l)(p+1) - 1]! ( 2 


2P+1 2(2"i-i)(p+i) 


am ^m—1 

n 

2 Clm—l ^m—1 

(2to - l)(p + 1) - 1 


p! [(2 to- l)(p+ 1) - 1]! 

22'"(P+1) [(to- l)!]2(p+ 1)2("S-1) ^ 


(C4) 

(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 


(CM) 
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where 


2(m- l)(p+ 1 )[to(p+ 1) - 1] 
^ (2p+l)[(2m-l)((p+l)-l] ■ 


(CIS) 


Note that if m » 1, 




m{p + 1) 
2p+l ’ 


or 




1 + 


(P+ 1) 

2(2p+ 1) 


Stirling’s formula, 


n! - n<l, 


(C16) 


(C17) 


(CIS) 


may be used to write 


{Am-By. - (^1 - 

/ Am \ "4"i-S+l/2 

e-«+i/2 


In the last step, we used. 


and 


We may show that 


B {Am)/B — 


^—{Am—B) 


lim 1 — 


1 - 


— ) 

Am J 


-B+1/2 


1 . 


"" 7- (P + l)2"^(P+l)-P-3/2 pP(l-2m) + l/2 + l)]p 

[[m — ij!J Att 


by application of Eq. (C19) to the factorials on both sides of this relation. 
This allows us to rewrite Eq. (C14) as 


(C19) 


(C20) 


(C21) 


(C22) 




f ^ 

, P+1/2 

r(p+i)pi 

VP + l/ 

) 

2PpP 


X [2pm — (p + 1)]!. 


(C23) 


For the case p = 3, this leads to Mn oc (3n — 4)!, which is the rate of growth found in Ref. [6]. This rate of growth in 
turn leads to the asymptotic behavior of the probability distribution given in Eq. ([^. 
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